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Abstract 1 

Two equation and higher order closures for compressible turbulence fail to capture the compressible 
wall layers’ log scaling. Accounting for the distinction between Favre and Reynolds averaged vari- 
ables in the compressible moment equations indicate that turbulent transport expressions obtained 
using the "variable density approximation 1 ’ are in error. The error is related to the enstrophy. a 
Reynolds averaged variable appearing in the equation for the Favre averaged k ; recognizing this fact 
an expression for the transport of dissipation consistent with simple mixing length arguments is 
obtained. Within the (limited) context of a gradient transport hypothesis a rational form for the 
turbulent transport of the dissipation is found. Modestly better agreement with the well established 
compressible Van Driest log scaling is found in a k — e calculation. 


] This research was supported by the National Aeronautics and Space Administration under NASA Contract No. 
NASI- 19480 while the author was in residence at the Institute for Computer Applications in Science and Engineering 
(ICASE), NASA Langley Research Center. Hampton, VA 23681-0001. 
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1. Introduction 

Huang et al. [1], Wilcox [2] have indicated that current k-e turbulence closures do not adequately 
capture the compressible log layer behavior. Part of the problem with these sorts of turbulence clo- 
sures appears to result from not adequately recognizing the distinction between Favre and Reynolds 
averaged quantities - in the Favre averaged equations both Favre and Reynolds averaged variables 
naturally appear. When mean density gradients are important, not accounting for the distinction 
contributes to poor results. 

In this article a derivation for the modelled dissipation equation for a compressible turbulent flow 
is described. While not accounting for many of the complex effects of compressibility, discussed 
in a general sense in Lele [4] or in the context of the fluctuating dilatation eg. Ristorcelli [5], the 
present development, which involves no modeling assumptions, shows better agreement with the 
compressible Van Dreist log law scaling. The usual gradient transfer expression for the turbulent 
transport of the dissipation is found to be written as utS/(pe s ) and not as n t pVe s as obtained using 
the so-called “variable density extension”. Huang [3] has communicated that such a form for the 
turbulent transport will produce better agreement with the Van Driest scaling. Wilcox [2] also 
comments on the fact that, if pe were the dependent variable, then k-e models would have closer 
agreement with the Van Driest scaling. While there has been speculation, Huang [1] about the 
properly conserved flow variable, there does not appear to be a clear indication as to what the 
proper formulation for the dissipation equation should be. This note indicates the appropriate 
form using a simple, but careful mathematical development. 

2. Derivation 

The nomenclature is now described. Upper case letters will be used to denote mean quantities and 
lower case letters fluctuating quantites. Exceptions to this rule are the mean density, p , (p has no 
convenient upper case form) or the mean viscosity p. Quantities with an asterisk denote the total 
field, mean and fluctuating: p* = p + p' or u* = U t + u t . Favre velocities will be denoted using 
the set [V,, u,] while Reynolds variables are denoted using the set [f7 t , u t ], The dependent variables 
are decomposed according to u’ = U; + = V] + r t - where < w, >= 0, {»,-} = 0, and p m = p + p' 

where < p' >= 0. The averaging operation is indicated using the angle brackets for time means. 
< ViVj >, and the curly brackets for the density- weighted or Favre mean {viVj}; the two averages 
are related by p{viV : } =< p*v t Vj >. Without loss of clarity the prime on the fluctuating density is 
dropped. The portion of the second-order moment equations, 

_ D , 

{ V z V j } — Pij 4- n ij + Tijk.k ~ < U j’}p (J ip > ~ < Ui,pCTj' p >, ( 1 ) 

of interest are associated with the viscous dissipation type terms, where = fi[u x ,j +Uj.i -2/3 u q . q 
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The quantities P tj , II, j, T, jk , k are, respectively, the production pressure strain and turbulent trans- 
port of the Favre averaged Reynolds stress, {vjVj}. Note that in the Reynolds stress equations the 
terms arising from the surface forces terms appear naturally in Reynolds, it,, variables while the 
problem is posed in Favre, v,-, variables. Repeated application of the product rule for differentiation 
and the definition of vorticity u>,- = e,j k Uj, k produces in the equation for the kinetic energy of the 
turbulence the usual dissipation quantities. Keeping only terms of interest produces 


P 


D 

Dt 


k = 


-p{viVj}Vi,j + T k - n < uj k u k > - -p < dd> 


( 2 ) 


where k = {vjVj}. The solenoidal dissipation has been rewritten in terms of the enstrophv assuming 
small scale isotropy. The dilatational dissipation is denoted e c = | < dd >. Our interest is with the 
portion of the dissipation associated with the vortical modes of the flow. It is conventional to define 
the solenoidal dissipation in terms of the Reynolds averaged enstrophy: pe s = p < u k u k >. In the 
“variable density approximation” and high Reynolds log layer limit the dissipation is assumed to 
be described by the high Reynolds number form of the Favre averaged equation 


Pjjj S s - [pvt^e 1 £s,q ],<? = P{CelPk ~ C e2 £*] £s/k. (3) 

where P k = S tJ = |[V' U +Vj„- -\Db x p in which D = U p , p is the mean dilatation which 

is small in the log layer portion of this flow. Note the location of the mean density in the turbulent 
transport term. The enstrophy’s replacement by what is treated as a Favre averaged variable, 
e s , must be done carefully. The conservation equation for the enstrophy, the primitive Reynolds 
averaged variable appearing in the Reynolds stress equations, is, to lowest order, 

< Up' + Uq Up + < UqUl >„j = — ^ ^ D - |- 2 < U>{S{jU>i !> jl < Uj , k , k > (4) 


The last two terms in the above equations are modeled by the right hands side of equation 3. The 
gradient transport hypothesis for the enstrophy transport produces < u q u 2 >= -v t < u P >,,, 
where u t - cjpjz. Replacing the enstrophy with the dissipation using p£ s = P < uj k u k > produces 

p—- t £ - [VtV7 l )„ ], g = ~\p eD + p[ c e\Pk - Ce2£s] £s/k. (5) 

The substantial derivative is along the Favre mean streamline. Note the location of the mean 
density in the turbulent transport term. In general the appearance of the mean density inside the 
first derivative will only be of importance when the mean density gradients are large, ^ ~ 

For weak density gradients the variable density assumption and the present derivation will give 
the same computational result. Calculations, to be shown shortly, were performed with a standard 
k - £ model with the usual gradient transport models for transport. The following values for the 
empirical constants were used: = 0.09, c e i = 1.44, c e2 = 1.92, at = 1.0, o z = 1.17. 


2 



The better performance of the k - u type of turbulence models comes from, in part, the fact that 
u — e/k is not a specific (per unit mass) variable and the Favre form of the convection operator is 
the correct expression for its turbulent transport. 

3. Conclusion 

One can follow Huang et al.'s [1] mean density gradient contribution analysis to asses the impact 
of these new forms of the transport terms in the k and 5 equations. The coefficients multiplying 
the density gradient terms in their equation (20) are notably smaller using the present formulation. 
Wilcox [2] has conducted a similar analysis in which he arrives at a similar conclusion. The 
implication of both these analyses is that the deviation from the law' of the wall will be smaller 
with these changes to the transport terms. This is substantiated in Figure 1 in which a Mach 8 
boundary layer calculation using a k — e turbulence model is shown. The top fine corresponds to 
“variable density approximation'': the lower fine(s) which are virtually indistinguishable correspond 
to the empirical compressible log law and an incompressible calculation. The middle fines are 
wfith the new r form of the transport in the dissipation equation given above: it show^s a modest 
improvement. The upper middle fine is without the bulk dilatation term. The point of this article 
is to present a rational form for the dissipation equation without adding ad hoc corrections. The 
compressible moment equations are, of course, complicated by several additional issues: the effects 
of the fluctuating dilatation and correlations involving the fluid property fluctuations, for example. 
These are the subject of current research and are expected to make additional contributions towards 
a better agreement with the w 7 ell established Van Driest scaling. 

To summarize: by careful accounting for the distinction between Favre and Reynolds averaged 
variables a different expression for the turbulent transport of the dissipation has been derived. The 
expression is consistent w r ith mixing length arguments adequate for w T all layers. Calculations for the 
Mach 8 bounday layer show that following a rational procedure produces a closer agreement with 
the log law r of the wall using the conventional models for the other terms in the modeled dissipation 
equation. 
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